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We present lattice QCD results on the neutron tensor charges including, for the first time, a 
simultaneous extrapolation in the lattice spacing, volume, and light quark masses to the physical 
point in the continuum limit. We find that the “disconnected” contribution is smaller than the sta¬ 
tistical error in the “connected” contribution. Our estimates in the MS scheme at 2 GeV, including 
all systematics, are g%T u = 1.020(76), g? = 0.774(66), g% = -0.233(28), and g 3 T = 0.008(9). The 
flavor diagonal charges determine the size of the neutron electric dipole moment (EDM) induced by 
quark EDMs that are generated in many new scenarios of CP-violation beyond the Standard Model 
(BSM). We use our results to derive model-independent bounds on the EDMs of light quarks and 
update the EDM phenomenology in split supersymmetry with gaugino mass unification, finding a 
stringent upper bound of d n < 4 x 10 -28 e cm for the neutron EDM in this scenario. 
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Low-energy precision measurements of neutron prop¬ 
erties provide unique probes of new physics at the TeV 
scale. Searches for the neutron permanent EDM d n 
have high sensitivity to new beyond the standard model 
(BSM) CP-violating interactions. Similarly, precision 
studies of correlations in neutron decay are sensitive to 
possible BSM scalar and tensor interactions. To fully 
realize the potential of the vibrant existing experimen¬ 
tal neutron physics program |Tj, one needs to accurately 
calculate matrix elements of appropriate low-energy ef¬ 
fective operators within neutron states. In this paper we 
describe lattice QCD calculations of the neutron tensor 
charges. In the future, these charges will be extracted 
with competitive precision from various measurements 
of the quark transversity distributions at JLab [2], and 
provide robust tests of the lattice results. 

The flavor diagonal charges g^’ d ’ s are needed to 
quantify the contribution of the quark EDM to the 
neutron EDM and thus set bounds on BSM sources 
of CP violation. We find that the contribution of the 
“disconnected” diagrams to g^, g j. and g f. are small. 
Our results on these charges allow us to constrain split 
supersymmetry models. 

The isovector charge g^~ u is needed in the analysis of 
precision neutron /3-decay. In Ref. 2] w e showed that to 
complement experimental measurements of the helicity 
flip contributions to neutron /3-decay at the precision of 
planned experiments (10 -3 level), we need to calculate 
the iso-vector scalar and tensor charges, g^~ u and g ^~ u , 
to about 10% accuracy. Results for g^T u presented here 
meet the desired accuracy with control over all systematic 
errors, while gg~ u requires 0(10) more statistics. 

Details of the lattice QCD calculations are given in 
a companion paper [4]. Here we summarize the main 
points and focus on the results using nine ensembles of 


Nf — 2 + 1 + 1 flavors of highly improved staggered 
quarks (HISQ) [5] generated by the MILC Collabora¬ 
tion j(I and described in Table [I] On these ensembles, 
we construct correlation functions using Wilson-clover 
fermions, as these preserve the continuum spin struc¬ 
ture. To reduce short-distance noise, all lattices were 
“HYP” smeared [7]. Extensive tests were carried out on 
these nine HYP smeared ensembles to look for the pres¬ 
ence of exceptional configurations [8], a possible problem 
with this mixed-action, clover-on-HISQ, approach. None 
were detected. Issues of statistics, excited state contam¬ 
ination, operator renormalization, lattice volume, lattice 
spacing and the chiral behavior are detailed in [3]. 

The flavor diagonal neutron charges g p are defined by 
(n(p, s)\0£\n(p, s)) = g^u s (p)Tu s (p) , with = qTq 
and the spinors satisfying J2 S w s(p)w(p) = (ft + m). 
The interpolating operator we use to create/annihilate 
the relativistically normalized neutron state | n(p, s)) is 
X (x) = e abc [qf(x)C 1 4(l+'y 4 )q b 2 (x)\ q{(x) with color 
indices {a, b, c}, charge conjugation matrix C, and c/i, c /2 
the two different flavors of light quark fields. 

The zero-momentum projection of xi x ) couples to the 
ground state, all radially excited states of the neutron, 
and multiparticle states. To reduce the coupling to 
radially excited states we Gaussian smear the quark 
fields in x( x )- To isolate the remaining excited state 
contamination, we include two states in the analysis 
of the two- and three-point functions at zero momen¬ 
tum [3]. Even though the excited state contribution is 
exponentially suppressed, we were able to isolate the 
leading two unwanted matrix elements (0|C?r|l) and 
(l|0r|l), where |0) and |1) represent the ground and first 
excited neutron states. We find that the magnitude of 
(O|0r|l) is about 16% of (0|C?r|0) and is determined with 
about 20% uncertainty on all the ensembles, whereas 
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Ensemble ID 

a (fm) 

Mr a (MeV) 

AC (MeV) 

L 3 xT 

M n L 

tsep / CL 

N conf 

N meas 

al2m310 

A 

0.1207(11) 

305.3(4) 

310.2(2.8) 

24 3 x 64 

4.55 

{8,9,10,11,12} 

1013 

8104 

al2m220S 

< 

0.1202(12) 

218.1(4) 

225.0(2.3) 

24 3 x 64 

3.29 

{8,10,12} 

1000 

24000 

al2m220 

0 

0.1184(10) 

216.9(2) 

227.9(1.9) 

32 3 x 64 

4.38 

{8,10,12} 

958 

7664 

al2m220L 

V 

0.1189(09) 

217.0(2) 

227.6(1.7) 

40 3 x 64 

5.49 

10 

1010 

8080 

a09m310 

<> 

0.0888(08) 

312.7(6) 

313.0(2.8) 

32 3 x 96 

4.51 

{10,12,14} 

881 

7048 

a09m220 

☆ 

0.0872(07) 

220.3(2) 

225.9(1.8) 

48 3 x 96 

4.79 

{10,12,14} 

890 

7120 

a09ml30 

o 

0.0871(06) 

128.2(1) 

138.1(1.0) 

64 3 x 96 

3.90 

{10,12,14} 

883 

7064 

a06m310 

□ 

0.0582(04) 

319.3(5) 

319.6(2.2) 

48 3 x 144 

4.52 

{16,20,22,24} 

1000 

8000 

a06m220 

o 

0.0578(04) 

229.2(4) 

235.2(1.7) 

64 3 x 144 

4.41 

{16,20,22,24} 

650 

2600 


TABLE I. The parameters of the (2+1+1) flavor HISQ lattices are quoted from Ref. [6j. The symbols used in the plots are defined 
along with the ensemble ID. All chiral analyses are carried out with respect to the clover valence pion masses M n which are tuned 
to be close to the Goldstone HISQ pion masses A/)) ea . We also give the source-sink separations (t se p/a) simulated, configuration 
analyzed ( N con {) and the total number of measurements (A^ meas ) made. Finite volume analysis is done in terms of M n L. 


|(l|C>r|l)| ~ (O|0r|O), but has 0(100%) errors. As illus¬ 
trated in Fig. [T]for the ci09m310 ensemble, the overlap of 
data in the center of the fit range for all the source-sink 
separations t sep indicates that excited state contamina¬ 
tion in the tensor charges is small and under control. 

The disconnected diagrams are estimated using a 
stochastic method accelerated with a combination of the 
truncated solver method (TSM) [5J [TO] . the hopping 
parameter expansion (HPE) [LQ jT5] and the all-mode¬ 
averaging (AMA) technique [13; . In most cases, the dis¬ 
connected contribution is small and consistent with zero 
as illustrated in Fig. |T] for the a09m310 ensemble. This 
feature was also observed in Ref. |Hj. We find that the 
light quark contribution is too noisy to extrapolate to the 
continuum limit, so we do not include it in the central 
value. We, however, use the largest estimate, 0.0121, on 
the coarsest ensemble al2m310 as an additional system¬ 
atic error in g j,, and g^, + “. 

The renormalization factor, calculated nonperturba- 
tively in the RI-sMOM scheme [Tol [TTii using the iso¬ 
vector operator, contributes a significant fraction of the 
total error. The charges converted into the MS scheme 
at 2GeV are given in Table [II] and Fig. [2j They are 
essentially flat in the three variables, lattice spacing a, 
the pion mass M n and the spatial lattice size L . We 
make a simultaneous fit to the data using the lowest or¬ 
der ansatz appropriate to our not fully O(a) improved 
clover-on-HISQ formulation: 

gria, M n , L) = ci + c 2 a + c 3 Af 2 + C 4 e~ MwL . (1) 

As discussed in [3], with current data the extrapolation 
to the physical point (M T = 135 MeV, a = 0, M n L = oo) 
is insensitive to additional corrections. The final renor¬ 
malized charges for the neutron E3 are 

g£ = 0.774(66), g £ =-0.233(28), 

g d~ u = 1.020(76), g^ +u = 0.541(67). (2) 

The y 2 /dof for the fits are 0.1, 1.6, 0.4 and 0.2, respec¬ 
tively, with dof = 5. Including the leading chiral loga¬ 


rithms [ 151 in Eq. 0 gives similar results [£ . g j,, after 
extrapolation in the lattice spacing a and Af 2 , is 

9t = 0.008(9), (3) 

with a y 2 /dof = 0.29 with dof = 2. The intercept of the 
fit on the [g|,, a] plane is shown in Fig. [3] 

Our result for g^ u , with control over all systematic 
errors, is in good agreement with other lattice calcula¬ 
tions ,11, 50]. The LHPC [51] and RQCD [55] Collabo¬ 
rations also find no significant dependence on the lattice 
spacing and volume, but do find a small dependence on 
the quark mass, so they extrapolate only in the quark 
mass using linear/quadratic (LHPC) and linear (RQCD) 
fits in AT 2 . Their final estimates, g^~ u = 1.038(11)(12) 
(LHPC) and g*T u = 1.005(17)(29) (RQCD) are consis¬ 
tent with ours. A fit to our data versus only Af 2 , shown 
as an overlay in Fig. [2] (center), gives a similarly accurate 
estimate g^T u = 1.059(29) with a y 2 /dof = 0.3. 

Our results on the tensor charges have implications 
for the neutron EDM and CP-violation in BSM theories. 
At the hadronic scale, g ~ 0(1) GeV, after integrating 
out all heavy degrees of freedom the dominant effect of 
new CP-violating couplings in BSM theories is encoded 
in local operators of dimension five and six. Leading, 
among them, are the elementary fermion EDMs [551 151] : 

SjC cpv =3 - | £ df fa^F^f . ( 4 ) 

f=u,d,s,e 

The contribution of the quark EDM d q to d n is [25, 26] 
d n = 9t du + 9t dd + 9t ds > (5) 

consequently, improved knowledge of g j, combined with 
experimental bounds on d n provides stringent constraints 
on new CP violation encoded in d q . 

Our calculation has the following impact: (i) We re¬ 
duce the uncertainty on g^ from the ~ 50% of previ¬ 
ous QCD sum rules (QCDSR) estimates [27: to the 10% 
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FIG. 1. Fits illustrating the excited state contribution in the connected g^~ u and light and strange disconnected diagram for 
the a09m310 ensemble. The data points represent <//(r, t sep ) obtained from calculations at different source-sink separations 
tsep and operator insertion times r. The solid black line and the gray band are the ground state estimate and error. 



a (fm) M n 2 (GeV 2 ) M n L 

FIG. 2. A simultaneous fit of neutron data versus a, M .and M n L using Eq. 0- The error band is shown as a function 
of each variable holding the other two at their physical value. The data are shown projected on to each of the three planes. 
The symbols are defined in Table [I] The extrapolated value is marked by a red star. The thin gray band and the dashed line 
within it in the middle panel show the fit versus M 2 assuming no dependence on the other two variables. 



FIG. 3. The data for gj and intercept of the fit versus a and 
M n on the [<?T,a] plane. Notation is the same as in Fig. [5] 



FIG. 4. Bounds on d U: d, defined in the MS scheme at 2 GeV, 
with a la slab prior on g£ and gtjn given in Eq. 0 and g J = 0. 


level. (For a comparison of the lattice results with the 
Dyson-Schwinger [28] and other methods [291432] see [4j.) 
(ii) The central values of g^' d are roughly 3/5 of the 
QCDSR and quark model estimates m widely used in 



FIG. 5. Iso-level lines of d n and d e in split SUSY in the Mi-g 
plane using sin^ = 1, tan/3 = 1, and central values of g^’ d,s . 

phenomenological studies of BSM CP violation, (iii) 
Bounding the strangeness tensor charge g^ at the per¬ 
cent level is important for a large class of models in which 
d q oc Trig since m s /md ~ 20. Our results imply that in 
such models g^d s may contribute up to 35% of the to¬ 
tal d n and the current 0(1) fractional uncertainty in gj, 
gives rise to the largest uncertainty in d n . The contri¬ 
bution of EDMs of heavier quarks to nEDM appears at 
two-loops and does not grow with m q . In this work, we 
ignore contributions of the charm (not calculated) and 
heavier quarks. 
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Ensemble ID 

con.d 

9t 

con.u 

9t 

con,d — u 

9t 

con ,d-\-u 

9t 

disc.Z 

9t 

disc.s 

9t 

al2m310 

0.852(37) 

-0.215(12) 

1.066(46) 

0.637(31) 

-0.0121(23) 

-0.0040(19) 

al2m220S 

0.857(43) 

-0.209(19) 

1.066(50) 

0.649(44) 

— 

— 

al2m220 

0.860(40) 

-0.215(15) 

1.075(48) 

0.644(36) 

-0.0037(40) 

-0.0010(27) 

al2m220L 

0.840(37) 

-0.194(12) 

1.033(45) 

0.647(33) 

— 

— 

a09m310 

0.840(28) 

-0.2051(98) 

1.045(34) 

0.634(25) 

-0.0050(22) 

-0.0005(21) 

a09m220 

0.836(28) 

-0.216(10) 

1.053(34) 

0.619(25) 

— 

-0.0021(54) 

a09ml30 

0.809(40) 

-0.222(20) 

1.032(44) 

0.587(45) 

— 

— 

a06m310 

0.815(29) 

-0.199(10) 

1.015(34) 

0.617(27) 

-0.0037(65) 

-0.0005(55) 

a06m220 

0.833(52) 

-0.264(22) 

1.099(59) 

0.569(55) 

— 

— 


TABLE II. Renormalized estimates of the connected (g“ n ) and disconnected (pt' sc ) 


contributions in the MS scheme at 2 GeV. 


While, in general, BSM theories generate additional 
CP-violating operators in Eq. ([4]), there exist models in 
which the fermion EDMs are the dominant sources of 
CP violation at low-energy, thus controlling the pattern 
of hadronic and atomic EDMs. For such cases, using 
Eq. (§, our results on the tensor charges, and the 
experimental limit on the neutron EDM (33] . we show 
90% confidence level (CL) bounds on quark EDMs d u ^ 
in Fig. [4] [511 . 

One notable scenario in which fermion EDM oper¬ 
ators provide the dominant BSM source of CP viola¬ 
tion is “split SUSY” (55H5T1 . in which all scalars, except 
for one Higgs doublet, are much heavier than the elec- 
troweak scale. This SUSY scenario achieves gauge cou¬ 
pling unification, has a dark matter candidate, and avoids 
the most stringent constraints associated with flavor and 
CP observables mediated by one-loop diagrams involving 
scalar particles. Contributions to fermion EDMs arise at 
two loops due to CP violating phases in the gaugino- 
Higgsino sector, while all other operators are highly sup¬ 
pressed [33 [33] ■ To illustrate the impact of improved 
estimates of matrix elements in split SUSY, we use the 
analytic results and setup of Ref. [38] . namely unified 
framework for gaugino masses at the GUT scale and a 
scalar mass m = 10 9 GeV. The light fermion EDMs 
d e ,u,d,s depend on a single phase </>, on tan/3 [approxi¬ 
mately through the overall factor sin(</>) sin(2/3)] and on 
the gaugino (M 2 ) and Higgsino (/r) mass parameters. 
Following Ref. [35] we set tan fi = 1, sin 0=1 and present 
the results as contours in the M 2 -/i plane, in the range 
between 200 GeV and 10 TeV. 

Figure [5] shows iso-level curves of d n as well as the 
curve d e = 8.7 x 1CU 29 e cm, corresponding to the cur¬ 
rent 90% C.L. limit [ID]. For the neutron EDM we use 
Eq. evaluating both the d g ’s and the tensor charges 
at the scale /r MS = 2 GeV. Our result for d n is apprecia¬ 
bly smaller (factor of ~ 3) than the one in Ref. [35]. We 
have traced back this difference to (i) our smaller values 
of the tensor charges compared to QCDSR n (ii) 
different values for the light quark masses: we use the 
PDG [HJ central value rrid( MS, g, = 2GeV) = 4.75 MeV, 
while the value corresponding to the quark condensate 



FIG. 6. Regions in \l 2 -g plane corresponding to d n /d e = 
2,3,4 in split SUSY, obtained by varying g^ ’ s within our 
estimated uncertainties. The lines correspond to d e = 8.7 x 
10~ 29 e cm for sin </> = 0.2,1. 

used in Ref. [ 31 ] is larger, rrid{ MS,/i = lGeV) ~ 9 MeV. 

In Refs. 33 S3, ^ was pointed out that the strong 
correlation between electron and neutron EDM would 
provide a valuable experimental test of split SUSY. To 
investigate this further, in Fig. [6] we present bands cor¬ 
responding to different values of d n /d e in the Mi~\i 
plane, the thickness of the bands reflects the behavior 
of d n /d e [M2 1 n] and the uncertainty induced by the ten¬ 
sor charges, dominated by g^. The fact that we can draw 
disconnected bands for d n /d e = 2,3,4 is a welcome con¬ 
sequence of our reduced uncertainties in using the 
QCDSR input, each band would be as thick as the whole 
plot, giving essentially no discrimination. 

Finally, based on the current 90% C.L. limit on d e , 
we derive an upper limit for the neutron EDM in split 
SUSY. By maximizing the ratio d n /d e along the iso-level 
curves d e = 8.7 x 10 -29 e cm corresponding to simp < 1, 
allowing g^, to vary in the lattice QCD ranges, we arrive 
at d n < 4 x 10~ 28 e cm [ 33 ] ■ Therefore, observation of 
the neutron EDM between the current limit of 3 x 10~ 26 e 
cm [ 33 ] and 4 x 10 -28 e cm would falsify the split-SUSY 
scenario with gaugino mass unification. 
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